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We present here a detailed granular flow characterization together with force measurements for the 
quasi-bidimensional situation of a horizontal cylinder penetrating vertically at a constant velocity 
in dry granular matter between two parallel glass walls. In the velocity range studied here, the drag 
force on the cylinder does not depend on the velocity Vb and is mainly proportional to the cylinder 
diameter d. While the force on the cylinder increases with its penetration depth, the granular 
velocity profile around the cylinder is found to be stationary with fluctuations around a mean value 
leading to the granular temperature profile. Both mean velocity profile and temperature profile 
exhibit strong localization near the cylinder. The mean fiow perturbation induced by the cylinder 
decreases exponentially away from the cylinder on a characteristic length A that is mainly governed 
by the cylinder diameter for a large enough cylinder/grain size ratio d/dg-. A ~ d/4 + 2dg. The 
granular temperature exhibits a constant plateau value To in a thin layer close to the cylinder of 
extension Stq '~ A/2 and decays exponentially far away with a characteristic length At of a few 
grain diameters (At ~ 3dg). The granular temperature plateau To that scales as V^dg/d is created 
by the fiow itself from the balance between the "granular heat" production by the shear rate Vb/A 
over Stq close to the cylinder and the granular dissipation far away. 

PACS numbers: 45.70.-n, 45.50.-j, 83.80.Fg 
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INTRODUCTION 



The characterization of forces on moving objects in 
granular matter is important in fields ranging from fluid 
mechanics to geophysics and biophysics vifith, for in- 
stance, the practical situations of meteoritic impacts on 
planets or asteroids [ij and of the motion of living organ- 
isms in sand [2] . As a matter of fact, a better understand- 
ing of meteoritic impacts requires a better knowledge of 
the forces experienced by impactors, and this gave rise to 
numerous physical studies at the laboratory scale 
and references herein) or by numerical simulations 
Most of these studies have identified mainly two terms 
in the forces experienced by an impactor: One term pro- 
portional to the square of the velocity and independent 
of the depth during the initial stage of the penetration 
(at high velocity) and one term proportional to the pen- 
etration depth and independent of the velocity at the 
end of the penetration (at low velocity) !3i, 4, q]. The 
stopping time has been shown to display a non intuitive 
behavior, with a smaller value for larger impact velocity 
0, 13, Hi • These behaviors for the penetration of the im- 
pactor may depend on the packing fraction of the grains 
[31 ; the packing fraction also determines the grain ejec- 
tion, which can give rise to a spectacular upwards vertical 
jet for a fluidized packi ng Q or to an opening corona 
for a dense packing [lol . Illl j. quite similar to what is ob- 
served for splashes in water of droplets or solid spheres 
p^ . Moreover, the motion of certain living organisms in 
grains, such as sand snakes, is interesting to understand 
in order to, for example, create artificial robotics able to 
move in grains [2|. In these last phenomena, the cou- 



pling between lift and drag forces is important [13|. All 
these physical phenomena require a better knowledge of 
the flow characterization around moving objects in grains 
and the associated forces. In this work, we shall focus on 
low velocity intrusions. 

Many studies have been performed on the drag force 
experienced by objects in relative vertical or horizontal 
motion with grains, but fewer studies have been done 
concerning grain flow characterization. In dense granu- 
lar matter, the drag force measured in velocity-controlled 
experiments has been found not to depend on the veloc- 
ity at low velocities, and to be proportional to the surface 
area of the object and roughly to its depth [13] . As in 
classical fluids, the corresponding drag force has been 
shown to depend on the exact shape of the object p^ . 
This force may depend on the packing fraction [l6j . on 
the vibration of the grains , or on possible dynamical 
air effects [l8!| . Amarouchene et al. reported a kinematic 
study of the flow around a thin disk placed in a Hele-Shaw 
cell and submitted to a vertical granular chute flow [T9| . 
The perturbed streamlines were observed to be located 
in a parabolic shape near the object with a triangular 
region of non flowing or slowly creeping grains, clearly 
due to the stabilizing friction force from the two close 
walls [13, mi . The measured velocity profile exhibited a 
non linear but monotonic variation away from the object. 
Chehata et al. performed the same kind of experiment for 
a cylinder in a much larger channel, and they measured 
at the same time the velocity field using a particle image 
velocimetry (PIV) technique and the drag force experi- 
enced by the cylinder in the grain chute flow [22|. They 
reported measurements of the flow vorticity and veloc- 
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ity fluctuations that showed shear locahzation near the 
cyhnder, and they noticed that the drag force is indepen- 
dent of the grain velocity and proportional to the cylin- 
der diameter. In a horizontal bidimensional experiment 
of vibrated disks at different packing fraction close to 
jamming, the measured velocity field around an intruder 
dragged at a constant force also exhibited shear localiza- 
tion [23i] . In a horizontal quasi bidimensional experiment 
in which a thin but large cylinder is pulled relative to a 
dense packing of millimetric alumina beads, Takehara et 
al. [2JI investigated motion at high velocities for which 
the force is shown to increase quadratically with velocity, 
and they observed high grain fiuctuation motions and ve- 
locity perturbation in a narrow zone around the cylinder 
that increases with the cylinder diameter. In the case of 
a dilute fiow, the drag force is also observed to be pro- 
portional to the square of the velocity, and proportional 
to the cylinder diameter and to the effective density of 
the fluid as in classical hydrodynamics with a drag coef- 
ficient that increases from about 1 at low Kn to about 
2-2.5 at high Kn [1^, H^, where the Knudsen number Kn 
compares the mean free path in the gas to the object size. 

In a previous work [23], we reported key results for 
both the experimental measurements and a hydrody- 
namic modeling of the canonical problem of the flow 
around a circular cylinder and the associated drag force 
in the granular case. In the present paper, we report 
the detailed experimental study of the mean grain veloc- 
ity profile and fluctuations around a horizontal cylinder 
penetrating vertically in granular matter together with 
drag force measurements. In Sec. |lTl we present the ex- 
perimental set-up. Experimental results on the profiles 
of both the mean velocity and the velocity fiuctuations 
are presented in Sec. IIIIl whereas Sec. IIVI concerns force 
measurements. The results are discussed in Sec.lVlbefore 
the conclusion. 



II. EXPERIMENTAL SETUP 

The experiments consist in the vertical penetration of a 
horizontal cylinder at a given velocity Vq into a packing of 
grains. The grains are rather monodisperse sieved glass 
beads of diameter dg ranging from 0.5 to 4 mm with a rel- 
ative dispersion Adg/dg of about 10% around the mean 
value, and density pg — 2.5 10'^ kg m~'^, contained in a 
rectangular box of length L — 0.2 m, height H — 0.1m 
and width b — 40 mm (Fig. [!}. The granular medium is 
prepared by gently stirring the grains with a thin rod, 
and the surface is then fiattened using a straightedge. 
We have checked that this preparation leads to repro- 
ducible results with only small variations. The solid vol- 
ume fraction is $ ~ 0.62 characteristic of a dense gran- 
ular packing, and the density of the granular medium is 
thus p = pg^ ~ 1.5 10'^ kg m~^. The cell containing the 
granular packing can be moved up along a vertical trans- 
lation guide at a given velocity Vq ranging up to 50 mm 
s~^ by a step by step motor (Mavilor BLS-55). A steel 




FIG. 1. (Color online) Sketch of the experimental setup. 



cylinder of length b and of diameter d ranging from 10 to 
40 mm is first maintained above the grains in between two 
vertical glass walls. The length b of the cylinder is man- 
ufactured about 0.2 mm smaller than the width b of the 
cell, and the cell is free to move along the axial direction 
of the cylinder. In this way, the cylinder has only mini- 
mal mechanical contact and thus solid friction with the 
walls during its vertical displacement, and no grain can 
fit between the cylinder and the walls and jam the sys- 
tem. As the cell moves upwards, the cylinder penetrates 
gradually at a constant velocity in the granular medium. 
The entire dynamics of the grains during the penetration 
is extracted by recording their motion at the front glass 
wall using a high-speed video camera that can take up to 
1000 images per second in the full resolution 1024 x 1024 
pixels. The grains are lighted from the front and a black 
curtain is put behind the cell so that the grains appear 
in a white on black background with a good contrast. 
The images are taken at a sampling rate / adjusted on 
the velocity Vq to have / > Vo/dg {e.g. f = 50 Hz for 
Vq = 10 mm s^^ and dg — 1 mm) so that the largest grain 
displacement between two successive images is smaller 
than one grain diameter. The successive images are then 
analyzed by Particle Image Velocimetry (PIV) software 
(Davis, LaVision) to get the velocity fields of the grains 
as shown in Fig. [5] The size of the final interrogation 
windows used in the correlation technique is set typically 
to one grain diameter to have the best resolution con- 
sidering the discrete nature of the medium. The spatial 
resolution of the obtained velocity field is thus one grain 
diameter. As the cylinder and video camera are fixed in 
the laboratory frame of reference and the grains move up, 
the velocity field of the grains is measured in the frame 
of reference of the cylinder. The grain velocity far from 
the cylinder thus corresponds to the undisturbed upward 
imposed velocity Vq, whereas the velocity is disturbed 
close to the cylinder. Considering the discrete nature 
of the granular medium and the finite size of the PIV 
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FIG. 2. (Color online) Left: Typical instantaneous velocity 
field obtained by PIV measurements for a cylinder of diameter 
d = 20 mm penetrating in a packing of glass beads of diameter 
dg = \ mm, at the velocity Vo = 50 mm s~^. Right: Sketch of 
the cylindrical coordinates and used notations. 



correlation window corresponding to one grain diameter, 
the PIV calculation is distorted at a distance from the 
cylinder surface smaller than half a grain diameter so 
that velocity measurements will be reported only outside 
this zone for radial distances r from the cylinder such as 
r > [d + dg) /2. With a force sensor (FGP Instrumenta- 
tion FN3030) of range 50 N, we also measure in the mean 
time the drag force experienced by the cylinder during its 
penetration in the granular packing at the sampling rate 
2 kHz. The cylinder is related to the force sensor by a ver- 
tical thin rod of diameter 3 mm and length 10 cm. With- 
out grains in the cell, we checked that the friction force of 
the cylinder on the glass walls is totally negligible. With 
grains in the cell, we made sure that no grain was fitting 
into the small spacing between the cylinder sides and the 
walls and was distorting the force measurement during 
the cylinder penetration. In the following we present ve- 
locity and force measurements that have been done for 
a cylinder/grain size ratio in the range 10 ^ d/dg ^ 80 
and for a cylinder velocity 0.5 ^ Vq ^ 50 mm s~^ remain- 
ing smaller than {gdg^/"^ ^ 70mm s~^. In this velocity 
range, the flow regime is expected to be quasi-static. The 
box length L was chosen large enough [L/d > 5) to not 
play any significant role in the measurements [28| . We 
also restrict ourselves to penetration depth values zi, not 
too close to (zf, < iJ — d) so that there is no significant 
role played by the bottom wall [2^, [2^ , and H will thus 
not be a relevant parameter for the present results. 



III. MEAN VELOCITY PROFILES AND 
VELOCITY FLUCTUATIONS 

A. Mean velocity 

The instantaneous velocity field obtained at each time 
t by the PIV technique is decomposed within the cylin- 
drical coordinates: 

v(r, e, t) = Vr{r, e, t) Gr + vg{r, e, t) ee, (1) 
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FIG. 3. (Color online) All 230 successive instantaneous (a) ra- 
dial and (b) azimuthal velocity profiles normalized by the pen- 
etration velocity Vb taken at the radial location r = 15 mm at 
regular intervals (At = 20 ms) during the penetration interval 
20 < 2i, < 70mm (d = 20mm, dg = 1mm, Vb = 10mm s~^). 

( ) Time average value and (■ • ■) corresponding standard 

deviation. ( — ) Cosine and sine fits of data by Eq. ((J} with 
Ar ~ 0.5 and Ae ~ 1.4. 



where r is the distance from the cylinder center and 9 
is the angle from the direction of motion (6* > anti- 
clockwise) (Fig. [2]) . We have checked that the granular 
fiow is bidimensional: No transverse fiow along the y di- 
rection perpendicular to the glass side wall can be seen 
either from the side (no grain appears at the glass wall or 
disappears from the glass wall during the cylinder pen- 
etration) or from the top (no transverse roll motion can 
be seen at the free surface). In addition, the results pre- 
sented in the following will be for — 7r/2 < 9 < tt/2, 
i.e. for the upstream grain flow, and we have checked 
that the downstream flow is very similar, with thus an 
upstream-downstream symmetrical flow field character- 
istic of low Reynolds numbers flows. Figure |3] shows the 
successive instantaneous (a) radial velocity component 
Vr{r, 9, t) and (b) azimuthal velocity component vg(r, 9, t) 
taken at a given radial location r (here r — d/2 = 5 mm), 
as a function of 9 and normalized by the penetration ve- 
locity Vq. These velocity profiles are quite noisy as they 
are not averaged at all but they do not show any sys- 
tematic time evolution during the penetration. This ob- 
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servation holds regardless of the radial position and the 
run. The velocity profile around the penetrating cylinder 
can thus be considered as stationary and is now time- 
averaged during each penetration run. The correspond- 
ing time average v(r, 0) = (v(r, 6',i)}f of the instanta- 
neous velocity profiles appears as a thick dashed line in 
Fig. [3] and the corresponding standard deviation is in- 
dicated by the surrounding two thick dotted lines. The 
radial velocity Vr{r^ 0) is maximum near the direction of 
motion (9 ~ 0) and vanishes around the equatorial plane 
{9 ~ ±7r/2), and the opposite is true for the azimuthal 
velocity vg{r, 9), which is maximal near the equator and 
vanishes near the axis of motion. This suggests the usual 
simple 9 dependence for the two cylindrical velocity com- 
ponents: 

v(r,e')/Vo = -Ar{r) cos 9 + Ae{r)iiin9 eg. (2) 

The fits of experimental and vg data by cosine and 
sine functions are shown as thick solid lines in Fig. [3] and 
are rather close to the averaged profiles. Doing such an 
analysis for all the radial locations r allows us to obtain 
the radial functions Ar{r) and Ag(r)^ which prescribes 
the radial variation of Vr and vg. Such functions are 
plotted in Fig. Ufa) and show a strong flow localization 
around the penetrating cylinder. As a matter of fact, 
Ar{r) and Ag{r) tend toward the value 1, corresponding 
to an undisturbed velocity field as soon as r > 40 mm and 
thus r > 2d. The radial function of the azimuthal veloc- 
ity Ag{r) exhibits an overshooting maximum {Ag > 1) 
as a consequence of the incompressibility for the present 
two-dimensional (2D) fiow. This maximum is located 
very near the cylinder at r ~ 17 mm < d and is fol- 
lowed by an inflection point. The radial function Ar{r) 
of the radial velocity increases from zero at the cylin- 
der surface up to the asymptotic value 1 at large r, with 
an exponential law as shown by the straight behavior 
of 1 — Ar{r) in the semi- logarithmic plot in the inset of 
Fig.HKa). With this flnding, assuming an incompressible 
flow with div v — {l/r)d{rvr)/dr+{l/r)dvg/d9 — 0, and 
taking into account the slip for the tangential velocity at 
the cylinder surface vg{r = d/2) 0, we adopt thus the 
following expressions for fltting Ar(r) and Ag{r) data: 



Ar{T) = 



r -d/2 + \s 



1 — exp — 



, , , r - d/2 + Xs - X 

Ag{r) = 1 + '— exp 

A 



- d/2 

r - d/2 
A 



(3a) 
,(3b) 



where A is the characteristic length of velocity radial vari- 
ation, and the length characterizes the velocity slip 
tangential to the cylinder surface. As a matter of fact, 
Ag{r = d/2) — As/A so that there is some slip for A^ ^ 0. 
The slipping length Is defined as {dvg/dr){r = d/2) = 
vg{r — d/2)/ls is related to As by Ig = As/(2 — As/A). 
As experimentally Ag{r — d/2) — As/A ~ 1, we have 
Is ~ As. Thus the length As corresponds roughly here 
to the slipping length Is- The curves corresponding to 
Eqs. ([5a|) are plotted as solid lines in Fig. Ufa) with the 
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FIG. 4. (a) Radial dependence of the radial and azimuthal 
velocity components (•) Ar and (o) Ae, for Vo — 10mm s~^, 
d — 20 mm and dg = 1mm. ( — ) Best fits by Eqs. (|3ap with 
A ~ 6.4mm and As ~ 9.7mm for Ar{r) data, and A ~ 7.6mm 
and As ~ 8.1mm for Ae(r) data. Inset: 1-Ar{r) data in log- 
lin plot, (b) Radial dependence of the stream function, B{r), 
obtained either from (•) Ar or (o) Ag data, and ( — ) best fit 
from Eq. ((S} with A ~ 7.1mm and As — 9.2 mm. The grey 
zone corresponds to the cylinder interior (r < 10 mm). The 
dashed lines stand for (a) A^^e = 1 and (b) B{r) = r which 
corresponds to an undisturbed velocity field. 



best-fitting values A ~ 6.4 mm and As — 9.7mm for ^, (r) 
data, and A ~ 7.6mm and As — 8.1mm for Ag{r) data. 
The fitting curves are close to the data, and the A and As 
values found independently from Aj-ir) and Ag{r) data 
are close, which means that the present dense granu- 
lar flow is nearly incompressible. Note that as As val- 
ues are close to half the cylinder diameter the prefactor 
(r — d/2 + Xs)/r in Eq. ((3a|) is close to one, thus explain- 
ing the exponential behavior observed for \ ~ A^ [inset 
of Fig. Ha)]. 

As the flow can be considered to be bidimensional in 
the (r, 9) plane, one can introduce the scalar stream func- 
tion il){r, 9) related to the velocity field by vg = —dip/dr 
and Vr ~ {1 / r)dil) / d9 . Considering the 9 dependence for 
the velocity field of Eq. ([2]) , the stream function can be 
written as 

'ip{r,9) ^ -Vo B{r)iiin9, (4) 
where the radial function B{r) is related to the velocity 
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FIG. 5. Variations of the characteristic length A and As of 
mean velocity profile with (a, d) the velocity Vb (for d = 
20mm and dg = f mm), (b, e) the cylinder diameter d (for 
Vb ~ 5 mm s'^ and dg = 1mm) and (c, f) the grain size 
dg (for Vb — 5 mm s~^ and d = 40 mm). The symbols are 
experimental data and the solid lines correspond to (a) A = 
7mm, (b) A = 3mm + 0.2d, (c) A = 10mm+ 1.5dg, (d) A^ = 
9 mm, (e) As = 0.4d, and (f) As = 18 mm. 



radial functions A^. (r) and Ag (r) hy A^ — B/r and Ag — 
dB/dr. The B{r) curves obtained either from Ar{r) or 
Ag{r) data of Fig. |4lja) are shown in Fig. lUJb). These 
curves are close, meaning again that the present granular 
flow is close to being truly incompressible. The average 
experimental B{r) data can be fitted by the following 
empirical function deduced from Eqs. (l3al) : 



B{r) - {r-d/2 + Xs) 



1 — exp 



r-d/2 
A 



(5) 



For large radial distance r away from the cylinder (r > 
2d), B(r) is close to linear in r, which corresponds to 
the undisturbed flow. The interesting domain is thus the 
one close to the cylinder (r < 2d) where B{r) deviates 
significantly from r. The fit of the mean B{r) data of 
Fig. [UJb) by the previous equation leads to A = 7.1mm 
and As = 9.2 mm, which correspond to intermediate val- 
ues between the A and A^ values deduced directly from 
Ar{r) and Ag{r) data. In the following, the presented 
A and As values have been obtained by the fit of B{r) 
experimental data. 

The A and As variations with the penetration velocity 
Vq, the cylinder diameter d, and the grain diameter dg 
are shown in Fig. [S] One can see that neither A nor As 
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FIG. 6. Temperature field T{x, z) around a cylinder of diam- 
eter d = 20 mm penetrating with the velocity Vb = 10 mm 
s^^ in a packing of grains of diameter dg — 1 mm. 



depend on the cylinder velocity Vq [Figs.[5Ua,d)] with the 
nearly constant values A = 7 ± 1 mm and As = 9 ± 1 mm. 
The As values appear to be proportional to the cylin- 
der diameter d with the linear fit As = 0.4d very close 
to the data [Fig. EJe)], and no significative dependence 
upon the grain diameter dg can be seen in Fig. [5Uf ) with 
the nearly constant value As — 18 ± 2 mm correspond- 
ing again to about 0.4d. The A values increase linearly 
with the cylinder diameter d with a slope around 0.2 and 
a significative non zero value around 3 mm correspond- 
ing here to 3dg that may be extrapolated for vanishing 
d [Fig. Eljb)]. The A variations with the grain diameter 
dg shown in Fig. [SJc) are small but may increase from 
about 10 to 15 mm with an extrapolated value d/4 for 
vanishing dg and a slope around 1.5. In conclusion, these 
two characteristic lengths A and As do not depend on 
the penetration velocity Vq and are mainly governed by 
the cylinder diameter d with a possible weak effect of the 
grain size dg on A when the cylinder diameter is not much 
larger than the grain size. 



B. Velocity fluctuations 

During each penetration run, we observe that the ve- 
locity field is stationary with a well defined time average 
value and some erratic time fiuctuations. In the preced- 
ing section, we focused on the time average velocity field 
v(r, 0), and we now present the time fluctuations whose 
amplitude can be quantified by the granular temperature 
defined as T{r,6) = ((v(r,6',t) - v(r,6'))^)t. A typical 
temperature field T(r, 6) around the cylinder is shown 
in Fig. [S] in gray scale with black for the lowest value 
and white for the highest value. A zone of high temper- 
ature can be seen very near the cylinder with a radial 
extension of a few millimeters around the cylinder sur- 
face. The precise radial dependence of the temperature 
extracted along the = streamwise direction in front 
of the cylinder is plotted in Fig. [7] for different cylinder 



6 




20 
r (mm) 

FIG. 7. Radial dependence of the temperature profile T{r) 
for 9 = 0, with d = 20 mm, dg = 1 mm and different velocities 
(V) Vo = 0.4 mm s~\ (•) Vc) = 2 mm s~\ (□) Vb = 10 mm 
s^^ and (a) Vb = 50 mm s~^. Guidelines for the eyes with ( — 

) the slope value At = 3.7 mm and ( ) the plateau value To 

and radial extension Stq = 4 mm close to the cylinder surface 
at r = d/2. 



velocities Vq. One can see in this semi logarithmic plot 
that all these different profiles have the same shape with 
about a plateau value denoted To in a zone of extension 
Sto close to the cylinder (0 < r — d/2 < Stq) followed 
by an exponential decrease T ~ exp(— t/At) with the 
characteristic length At at larger r (r — d/2 > St,,)- The 
different curves at different velocities Vq show that Tq in- 
creases with Vo whereas neither Stq nor At varies signifi- 
cantly as the decreasing parts of T{r) are about parallel: 
6to — 4 mm and At — 3.7 mm. 

The variations of Tq and At with the velocity Vq, the 
cylinder diameter d, and the grain diameter dg are shown 
in Fig. m The log-log plot of Tq with Vb in Fig. M^a) ex- 
hibits a clear data increase with a slope 2 leading to the 
natural scaling Tq ~ Vq . Figures ID^b) and show 
that the typical temperature Tq decreases with increas- 
ing cylinder size d and increases with the grain size dg. 
The Tq dependence upon dg is linear and the Tq depen- 
dence upon d may be hyperbolic such that Tq may finally 



scale as VQ{dg/d), which will be supported by the plot of 
Fig.lT^a). as we shall see in the discussion (Sec. V). Con- 
cerning the characteristic length At of the exponential 
temperature decrease away from the cylinder, Fig. |Sl^d) 
shows that At does not depend on Vq with the nearly con- 
stant value At — (3.7 ±0.5) mm. One can see in Fig.[Sl[e) 
that At does not depend significantly on the cylinder di- 
ameter d with a value around At — 4 mm corresponding 
thus to a few grain diameters. This is corroborated by 
Fig.lSJf), where At appears to increase linearly with the 
grain diameter dg if one ignores the singular point at the 
smallest grain size dg = 0.5 mm. The dependence of the 
radial width of the temperature plateau 6t„ on Vb, d, 
and dg is not shown independently, but Fig. IH] displays a 
clear correlation of Stq with the characteristic length A of 
mean velocity variations. This means thus that Stq does 
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FIG. 8. Temperature plateau To and characteristic length At 
of radial temperature variation as a function of (a, d) the ve- 
locity Vo (for d = 20 mm and dg = 1mm), (b, e) the cylinder 
diameter d (for Vo = 5 mm s^^ and dg = 1 mm) and (c, f) the 
grain size dg (for Vo = 5 mm s~^ and d = 40 mm). The sym- 
bols are experimental data and the solid lines correspond to 
(a) To = 0.05V^, (b) To oc d'^ (c) To oc dg, (d) At = 3.7mm, 
(e) At = 3.5 mm-|-0.03ci, and (f) At = 3.3 mm -|-l.ldg. 
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FIG. 9. Radial extension Stq of the temperature plateau as 
a function of the characteristic length A of the mean velocity 
profile (same symbols as in Fig. 



not depend on Vq and is mainly governed by the cylinder 
diameter d with only a weak effect of the grain size dg in 
contrast with the At size dependence discussed above. 
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IV. FORCE ON THE CYLINDER 

The force on an object in relative motion with a gran- 
ular material has been extensively studied these last few 
years by many authors in a lot of geometries. The drag 
force is observed to vary linearly with the surface area 
of the object and increases with the depth z, F ~ z", 
with possible non linear effects (a ^ 1) depending on the 
direction of motion and geometry (14, 15, 22, 29 31]. For 
a vertical rod in a horizontal granular flow, a is signifi- 
cantly larger than 1 (a ~ 2 [l^,!!^]); which is linked to the 
linear increase of the immersed surface with z. For ob- 
jects whose immersed area does not vary with the depth, 
the exponent a is closer to 1. For the vertical penetration 
of spheres, cubes, and horizontal rods, a ~ 1.2 was found 
[3^ . For the vertical penetration of spheres and also of 
cylinders and cones with aspect ratio one and vertical 
axis, Ref. (3l| found that a. > l,with a a value depend- 
ing on the geometry: a ~ 1 for these cylinders, a ~ 1.3 
for spheres, and a ~ 1.8 for these cones. This supralin- 
ear depth dependence of the force (a > 1) is observed 
at small depth but is followed by a sublinear dependence 
(a < 1) at large depth attributed to the counter force 
generated by the filled-in grains on the top of the ob- 
ject For a horizontal plate penetrating a granular 
medium, Ref. [1^ found a linear increase of the force with 
depth (a ~ 1) at small depth, followed by a saturation at 
larger depth (a ~ 0) attributed to the Janssen screening 
wall effect, and finished by an exponential increase very 
close to the bottom wall. 



We present here force measurements associated with 
the velocity profiles presented above and corresponding 
to our geometry of a horizontal penetrating cylinder in 
vertical displacement in the granular packing. A typical 
force variation as a function of the penetration depth 
z\y is shown in Fig. llOf a). where is the position of 
the cylinder bottom, with thus = when = 0. 
Figure [TUf a) shows that F increases with Zh whereas the 
grain velocity profile around the cylinder has been shown 
to be stationary in section UlI Al The F{zh) variation is 
about linear, and if one considers a hydrostatic equivalent 
"pressure" p = pgz that varies linearly with depth in the 
granular material, this would lead to the Archimedean 
force Fa = pgnhcP/A equal to about 0.2 N for the case 
of Fig. [rUTa). which is thus negligible when compared to 
the measured drag force. 

We do not observe any drastic change in the rate of in- 
crease of F with Zf), which may be due to the fact that we 
do not explore very large depth values nor depth values 
close to the bottom wall. The mean slope of F{zb) can 
be extracted and is displayed in Fig. [10] for different (b) 
penetration velocities Vqj (c) cylinder diameters d, and 
(d) grain diameters dg. One can see that the drag force 
is approximately independent of Vq and proportional to 
the cylinder diameter d. The F variation with the grain 
diameter dg is more complex: When F seems to be about 
constant for large grain diameters [dg > 1mm), it seems 
to increase with decreasing dg for sma ll g rain diameters 
{dg < 1 mm) as already mentioned by [ij, |23] . 

V. DISCUSSION 
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FIG. 10. (a) Drag force F on the cylinder as a function of its 
penetration depth Zb for Vo = 10 mm s~^, d = 20 mm, and 
dg = 1mm. Depth force variation {F/zh), with (■) standing 
for an average of z^ in the range d/2 < zt < H — d, as a 
function of (b) the velocity Vo (for d = 20 mm and dg — 
1mm), (c) the cylinder diameter d (for Vo — 5mm s~^ and 
dg — 1 mm), and (d) the grain size dg (for Vo = 5 mm s~^ and 
d = 40 mm). Symbols are experimental data and the solid 
lines correspond to the best fits of the data with (b) (F/zt) — 
0.29 N mm"^ and (c) ~ 0.012 N mm"^ 



The velocity field and the drag force around a mov- 
ing cylinder in granular matter show important differ- 
ences compared to the case of a Newtonian fluid. As 
a matter of fact, the drag force in granular matter is 
independent of velocity in the explored range of veloci- 
ties but increases with depth. This differs from Newto- 
nian fluids for which the hydrodynamic drag increases 
with velocity but is independent of depth, as long as 
fluid viscosity is independent of pressure so that pres- 
sure only leads to a constant Archimedean buoyancy 
term. Concerning the velocity field, the perturbation 




(b) 

60 80 100 



FIG. 11. Dimensionless characteristic lengths of the velocity 
profiles, (a) X/dg and (b) Xa/dg, as a function of the size ratio 
d/dg. ( — ) Best fit of the data (same symbols as in Fig. [5]) 
with equations (a) A = 2.3dg -|- 0.24d and (b) As = 0.44d. 
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FIG. 12. Temperature plateau To as a function of V^dg/d. 
( — ) Fit of equation To — I.^Vq dg / d through the data (same 
symbols as in Fig. [S]). 



created by the cylinder in granular matter appears lo- 
calized near the cylinder with an exponential radial de- 
crease {B{r) ^ r[l — cxp(— r/A)]), which differs signifi- 
cantly from the Newtonian case in either the inviscid case 
for which the perturbation relaxes with a power law of the 
distance {B{r) = r(l — cP /Ar'^)) or the viscous case for 
which the perturbation relaxes slowly at large distance 
from the cylinder {B{r) ~ r ln((i/2r)/ In(Re) for Re <C 1 
and r » d). The present granular case is closer to yield 
stress fluids where a strong localization is also observed 
[I^IIJ]. All our A measurements at different velocities, 
and cylinder and grain diameters, collapse in Fig. IllT a) 
when made dimensionless with the grain size dg leading 
to a linear increase oi X/dg with the size ratio d/dg, with 
a small but non zero value 2.3 ± 0.5 for vanishing d/dg 
and slope 0.24 ±0.03. This means that the characteristic 
length A of mean velocity radial variation is mainly gov- 
erned by the cylinder diameter d as long as the size ratio 
d/dg is greater than about 10, which is always the case 
here in contrast to small intruder size situations studied, 
e.g., in In contrast. Fig. fTlTb) shows undoubtedly 

that the "slipping length" As depends only on the cylin- 
der diameter d with the scaling As ~ (0.44 ± G.4)d very 
close to d/2. 

As the typical granular temperature Tg of the granular 
flow around the cylinder is connected with the relative 
velocity Vq, there is clearly a strong coupling between 
the mean velocity profile and the temperature profile. In 
that sense, the present situation resembles the case of the 
motion of a hot cylinder or sphere in a fluid of tempera- 
ture dependent viscosity. This last situation has already 
been studied as it corresponds to the geophysical situa- 
tion of the ascending motion of hot diapirs [35, 361 ■ 
this last case, the flow is also strongly localized close to 
the hot sphere where the viscosity is much lower than far 
away, as heat diffuses into a fluid layer near the object, 
thus lowering the fluid viscosity. Indeed, stress continuity 
causes regions of low viscosity to be regions of large strain 
rates, producing a temperature-induced shear band. In 
this last case, heat is produced outside the flowing mate- 
rial by the hot sphere. In the present granular case, the 
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FIG. 13. Dimensionless characteristic length (a) of tempera- 
ture plateau Sr^/dg and (b) of temperature exponential de- 
crease Xr/dg as a function of the size ratio d/dg. ( — ) Fit of 
equation (a) Stq = l-4dg -I- 0.12d and (b) At = 2dg + 0.07d 
through the data (same symbols as in Figs. [8] and [9]). 



"heat" production is due to the flow itself, and the scal- 
ing of the temperature plateau Tq appears very close to 
To = {1.5±0.5)V^dg/d in the log-log plot of Fig.[12]where 
all the data for different Vq, d, and dg collapse along a 
slope 1. We have seen in the preceding section (Fig. ^ 
that the radial extension Stq of the temperature plateau 
close to the cylinder is linked to the characteristic length 
A of the radial variations of the mean velocity. The plot 
of all data of Sxo/dg at different Vq, d and dg as a func- 
tion of d/dg in Fig. llSf a). exhibits a small but non zero 
value 1.4 ± 0.2 for vanishing d/dg and a slope 0.12 ± 0.02 
corresponding roughly to Stq ~ A/2. As A, Stq is mainly 
governed by the cylinder diameter d as long as the size 
ratio d/dg is larger than about 10. All these findings 
mean that the granular heat is created by the shear flow 
Vq/X ^ Vo/d close to the cylinder, in a region of radial ex- 
tension ^ A (more precisely Stq ^ -^/2)- This shear flow 
leads to velocity fluctuations Vodg/X ^ V^dg/d at the 
grain scale that in turn produce the granular tempera- 
ture {Vodg/d)'^. The observed scaling for the temperature 
plateau Tq V/^dg/d should come from an equilibrium 
between the corresponding "heat" production by the high 
shear close to the cylinder and the granular dissipation 
away from the cylinder. The characteristic length At of 
the exponential decrease of the temperature observed far 
away from the cylinder is mainly governed by the grain 
size dg-. all the data of Xr/dg at different Vq, d and dg 
plotted in Fig. I13f b) as a function of the size ratio d/dg 
only show a slight increase with d/dg of slope 0.07 ±0.07 
from the value 2 ± 1 at vanishing d/dg = 0. 

The flow localization, observed in our case in the vicin- 
ity of the cylinder, is indeed reported in various physical 
systems. It is a common feature of matter with granu- 
larity as recently reviewed by [st']. Such shear bands are 
currently observed for flow of granular materials, foams, 
or emulsions. In such systems, the thicknesses of the 
shear bands are usually of few grains or bubbles. Despite 
the generality of such observations, there is still no clear 
unifying framework for such shear bands 37[ . In the case 
of yield stress fluids, the shear-bands widths may be large 
compared to the microscopic sizes (38j . and shear local- 
ization is observed in Newtonian fluid flow close to walls 
at temperatures larger than the temperature of the fluid 
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far away (35l.[36j. 

Note that in a cylindric Couette device, Refs. [3^ \^ 
have reported a velocity profile for the sheared dry grains 
independent of the cylinder angular velocity and also in- 
dependent of the "pressure" imposed by a possible exter- 
nal upward or downward air flow. This has to be related 
to our case in which the velocity profile around the mov- 
ing cylinder does not depend on the depth and thus on 
the effective "pressure" . References [39, 40] also reported 
a small zone of constant granular temperature Tq ~ Vq 
and extension 6to — 2.8dg close to the rotating cylinder 
followed by an exponential decrease at larger distance 
with a characteristic length Ay — 4.7dg. The temper- 
ature profile we find in the present study for the fiow 
around a cylinder is similar, with a similar scaling of At 
but different scalings for Tq and Stq that should come 
from geometric consideration. 

In [2^ . a model was briefiy presented for the granu- 
lar flow around the cylinder based on continuous con- 
servation equations for mass, momentum, and granular 
temperature. The phenomenological equations relating 
normal and shear stresses, density, and temperature are 
given by the kinetic theory of granular gases [4l| in the 
dense limit. In this approach, the granular material may 
be viewed as a fluid whose viscosity decreases, at a given 
pressure, with the granular temperature. From a thermal 
point of view, the heat is produced by shearing the fluid, 
then diffuses into the material, and is dissipated during 
particles collisions. It follows from such an approach that 
a sheared zone increases the granular temperature, which 
fluidizes the material, increasing the ability to flow. The 
numerical simulation of mechanical and heat equations 
shows indeed the presence of such a fluidified zone near 
the obstacle ;27|. The constitutive equation relating the 
viscosity to the temperature and the pressure leads to 
a velocity-independent drag force in the creeping flow 
regime, as observed experimentally. 

Considering now our observations for the resulting 
drag force on the cylinder, one can infer an effective fric- 
tion coefficient /i = F/ (pgzbTrdb) based on a hypothetical 
effective hydrostatic pressure pgz at the depth z inside 
the granular material. All the data from different Vq, d, 
and dg collapse when plotted as a function of the size 
ratio d/dg in Fig. 1141 meaning that no other finite size 
effect (such as with the container sizes H, L, or b) is ob- 
served in the present experiment. The effective friction 
coefficient appears to be greater than 5 with a plateau 
value of about 7±2 at low enough d/dg {d/dg < 50) and 
a possible increase for larger d/dg values {d/dg > 50). 
This unusual value of the friction coefficient, which is 
about 20 times higher than the usual values of friction 
coefficients in dense granular fiows , has already been 
mentioned in [3| for the depth-dependent term of the drag 
force extracted from impact experiments. This means 
that the "pressure" is undoubtedly far from hydrostatic 
in the present situation even if the force variation with 
depth is almost linear. 
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FIG. 14. Normalized drag force F/ {pgztndb) on the cylin- 
der as a function of the size ratio d/dg. Same symbols as in 
Fig. [To] 



VI. CONCLUSIONS 

In this paper, we have studied both velocity profile and 
force measurements for the granular flow around a hor- 
izontal cylinder in the vertical penetration case. While 
the force increases with the depth, the velocity profile 
is shown to be stationary during the penetration with a 
well-defined average and some erratic fluctuations. The 
granular flow being close to bidimensional and incom- 
pressible, the stream function has been extracted from 
the velocity measurements. The flow perturbation is 
strongly localized close to the cylinder, exhibiting an ex- 
ponential radial decrease away from the cylinder with a 
characteristic length A that scales mainly with the cylin- 
der diameter d for large enough cylinder to grain size 
ratio d/dg-. X 2dg + d/4 « d/4 for d/dg > 10. The 
velocity fiuctuations quantified by the so-called granu- 
lar temperature T show also a strong localization near 
the penetrating cylinder, with a plateau value Tq in a 
narrow crown of extension Stq around the cylinder fol- 
lowed by an exponential radial decrease with a small 
characteristic length At- The scaling of Tq appears to 
be simply Tq ~ Vgdg/d coming from a balance between 
the "granular heat" production by the shear flow Vq/X 
Vo/d for d/dg > 10) over the distance Sto ^ A/2 
(~ d/8 for d/dg > 10) and the granular dissipation far 
away (r > Stq)- This granular dissipation which is at 
the grain scale leads to a characteristic length scaling At 
for the temperature decrease of a few grain diameters: 
At ^ 3c?g. 

In conclusion, the granular flow around a moving cylin- 
der is very different from the flow of a Newtonian fluid 
either in a viscous or inertial regime: For the Newto- 
nian fluid case, whatever the regime is, the drag force 
increases with the velocity and does not depend on the 
depth, whereas the contrary is observed for the granu- 
lar case for which the drag force does not depend on the 
velocity but increases with the depth. In an upcoming 
paper, we will detail a possible hydrodynamic modeling 
of this non classical fluid behavior based on kinetic theory 
adapted for dense dissipative granular systems as briefly 
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presented in [27[ . The shear locahzation characterized in 
the present experiment is certainly important to mider- 
stand the influence of a close wall on the motion of objects 
in some practical situations such as the dynamics of im- 
pacting spheres [28„ 43] . In addition, other experiments 
have to be performed to understand the very different 
drag force experienced by an object when moving down 



or up in penetration or withdrawal situations, as already 
observed by [s^l and [l6j . 
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